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Abstract 

We present an analysis of extensive large-scale Monte Carlo sim- 
ulations of self-avoiding fixed- connectivity membranes for sizes (num- 
ber of faces) ranging from 512 to 17672 (triangular) plaquettes. Self- 
avoidance is implemented via impenetrable plaquettes. We simulate 
the impenetrable plaquette model in both three and four bulk di- 
mensions. In both cases we find the membrane to be flat for all 
temperatures: the size exponent in three dimensions is v = 0.95(5) 
(Hausdorff dimension dn = 2.1(1)). The single flat phase appears, 
furthermore, to be equivalent to the large bending rigidity phase of 
non- self- avoiding fixed-connectivity membranes - the roughness expo- 
nent in three dimensions is £ = 0.63(4). This suggests that there is a 
unique universality class for flat fixed-connectivity membranes with- 
out attractive interactions. Finally we address some theoretical and 
experimental implications of our work. 
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1 Introduction 



There are many different physical systems occurring in nature that we can 
theoretically model as fixed- connectivity membranes |], ^, f| §]. A fixed- 
connectivity membrane is a two-dimensional mesh, in a d- dimensional bulk 
space, whose nodes have a time-independent coordination number. In other 
words the connectivity of the mesh is fixed in time, although there may be 
some spatial variation, provided the resultant density of defects (dislocations 
and disclinations) is small enough to preserve, at zero temperature, the quasi- 
long-range translational order and true long-range bond orientational order 
characteristic of two-dimensional crystals [0, || [7|. The statistical mechanics 
of such membranes is a two-dimensional version of the already rich statistical 
mechanics of essentially one- dimensional polymers |8], [], [3], |j. The physical 
bulk dimensions are d = 2, corresponding to monolayers or planar crystals, 
and d = 3, corresponding to membranes. The intrinsic importance of this 
universality class of membranes is enhanced by the wealth of examples pro- 
vided by the natural world. These include the spectrin/actin cytoskeleton of 
mammalian red blood cells ||10|| , polymerized amphiphilic bilayers and mono- 
layers fill , polymerized polymer sheets |12| , graphitic oxides [0, [14], [15| , 



metal dichalcogenides |TjJ and large sheet molecules of glassy B 2 Os |T7 
Experimental measurements of the physical properties of these real-world 
fixed-connectivity membranes can be compared with the results of theoret- 
ical models and numerical simulations. This is an active area of research 
on the experimental side, with the exciting possibility of novel membranes 
such as two-dimensional cross-linked DNA networks. A systematic and deep 
understanding of this field is also essential in developing theoretical models 
of more complex systems, such as the complete cell membrane, in which a 
cytoskeleton (a fixed-connectivity membrane) is coupled to an amphiphilic 
bilayer (itself a liquid-like membrane) by protein junctional complexes ^M . 

The statistical mechanics of phantom (self-intersecting) fixed-connectivity 
membranes is governed by the competition between the elastic energy of 
stretching and shearing deformations and the bending energy of shape fluc- 
tuations. This competition leads to a scale-dependent running of the bend- 
ing rigidity and elastic moduli of the Hamiltonian in the presence of ther- 
mal fluctuations, with the membrane stiffening under bending and softening 
under elastic deformations as the length scale grows. The upshot is a low- 
temperature ordered phase, characterized by an infinite persistence length 
and a spontaneous orientation of the membrane in the bulk space, and a high- 
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temperature crumpled phase with a finite persistence length and no preferred 
orientation in the bulk. The remarkable ordered phase, conventionally known 
by the inapt name flat phase, is one of the key reasons for the excitement 
in this field. Furthermore the global flat and crumpled phases are separated 
by a continuous crumpling transition associated with an ultraviolet-stable 
fixed point in the renormalization group flow of the bending rigidity. The 
crumpling transition is another exciting feature of the physics of phantom 
fixed-connectivity membranes. 

Realistic fixed-connectivity membranes are, however, self-avoiding (SA) 
in the sense that self-intersections will be energetically costly. This self- 
avoiding, or excluded-volume, interaction arises at the microscopic (nanome- 
ter) scale from the short range repulsion between any two monomers in the 
membrane. Study of this self-avoiding interaction has proven remarkably 
difficult. Analytical studies (see |3| and || for recent reviews) are ambigu- 
ous in that there is a variety of results for, say, critical exponents. This has 
made it difficult to know the precise nature of the global phase diagram and 
to reliably predict the key critical exponents. Numerical results using the 
balls and springs (BS) model, to be discussed in more detail below, clearly 
indicate that any degree of self-avoidance flattens the membrane at all tem- 
peratures. In other words the crumpled phase is unstable to the inclusion 
of self-avoidance. With its loss goes also the crumpling transition. It has 
been argued, however, that the loss of the crumpled phase is a result of the 
particular discretization of self- avoidance within the BS models, rather than 



a universal feature of self-avoidance itself |L9| . In particular distant-neighbor 
interactions due to the excluded volume of the balls can induce an effective 
bending rigidity for these models even with vanishing bare bending rigidity 
jl9| . Subsequent careful studies of BS models support the claim that the 
crumpled phase is lost (see Sec.0). More light can be cast on the problem, 
however, by exploring an alternative discretization of a SA fixed-connectivity 
membrane known as the impenetrable plaquette (IP) model. This will be the 
focus of the work presented in the rest of the paper. Impenetrable plaquette 
models have the advantage that the described membrane is extremely flexible 
in that individual plaquettes are completely free to fold on themselves. As 
a result they do not suffer from induced bending rigidity. The organization 
of the paper is as follows. In the next section we provide a quick overview 
of the BS models for completeness. In Sect.|| we introduce and discuss the 
impenetrable plaquette model studied in this paper. In Sect{| we describe 
the details of the algorithms used to numerically simulate the model. A 
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complete presentation of our results is given in Sect.[| for the cases of bulk 
dimensions d = 3 and d = 4. Finally we summarize in Sect.|6| and discuss the 
implications of our work for the experimentalist and outline some promising 
future directions. 



2 Self- avoidance: Numerical Studies 

The self-avoidance interaction at long distances (typically of the order of a 
micron) may be modeled by generalizing the Edwards model of a self-avoiding 
polymer [[2(J to objects with internal dimension D: 

UsA = \\ ^^VW - f(x')) , (1) 

where b is the dimensionless parameter governing the strength of self-avoidance. 

This implementation of self-avoidance captures the underlying universal 
physics for a broad class of microscopic potentials describing the interactions 
between the monomers of the membrane at the nanometer scale. 

The most direct discretization of fixed-connectivity membranes including 
self-avoidance is a network of iV monomers arranged in a triangular array. 
Nearest neighbors interact with a potential 

w / for |f| < p 

V NN {r) = < ', ' r , (2) 

I oo tor |r| > p v ' 

or a smoother variation. The tether length p is of the order of a few lattice 
spacings. The self-avoidance, or excluded volume, is introduced as a repulsive 
hard sphere potential, now acting between any two monomers in the mem- 
brane, instead of only nearest neighbors. A typical hard sphere repulsive 
potential is 

/ oo for | f| < a , s 

V **°W = { for \f\ > a ' (3) 

where a is the range of the potential, and a < p. Once again some smoother 
versions of this potential, continuous at |f| = a, have also been adopted. 
This model may be pictured as a network of springs, defined by the nearest- 
neighbor potential Eq.Q, with self-avoidance enforced by the finite radius a 
of the balls (EqJJ) and we will often refer to it as the BS (balls and springs) 
model. 
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There is conclusive evidence that the phase diagram of the three-dimensional 
BS model consists of a single flat phase - the crumpled phase being lost en- 
tirely [19,22-28]. A detailed review may be found in [[!]. 

Simulations for bulk dimension above three have also been performed 
22], |23| and provide clear evidence that the membrane remains flat in three 



and four dimensions and undergoes a crumpling transition for bulk dimension 
d > 5, implying that the lowest dimension, d c , in which a crumpled phase 
exists is five (d c = 5). 

It has been argued [|19| that the absence of a crumpled phase is due to 
induced bending rigidity from distant (next-to-nearest neighbor) excluded 
volume effects. This may place the self-avoiding membrane in the flat phase 
of the phantom membrane phase diagram (see and references therein), 
in which case one expects the flat phase of the self-avoiding and phantom 
membranes to be equivalent. If this is indeed the case, sufficiently weak 
excluded volume interactions should induce a bending rigidity below the 
critical crumpling transition coupling, leading to a crumpled self-avoiding 
phase. No conclusive evidence for this self- avoiding crumpled phase has been 
reported. This suggests that flatness is an inevitable consequence of self- 
avoidance combined with the fixed-connectivity (integrity) of the membrane. 
This picture would be much more convincing if it could be shown to apply 
to other models of self-avoiding fixed-connectivity membranes. 

An alternative model of self-avoidance is provided by the impenetrable 
plaquette model (IP), first simulated in P5| . These authors found a 
size exponent corresponding to a crumpled phase and compatible with the 
standard Flory estimate. A subsequent simulation |26j of a variation of the 
IP model, contradicted the results of [|24], ^] by claiming a flat phase, but 
found a fractal size exponent, from the finite-size scaling of the radius of 
gyration, of v = 0.87 (no errors quoted). This corresponds to a Hausdorff 
dimension dn = 2.3, which is neither Flory nor flat. In contrast, the analysis 



of the orientationally-averaged structure function in |26] gives a different 
fractal value v = 0.75(dn = 2.67) (no errors quoted). In short the status of 
the phase diagram for the IP model is currently murky, in contrast to the 
clear picture emerging from the BS models outlined earlier. 

Motivated by this unsatisfactory state of affairs we turn now to an analysis 
of our particular realization of the IP model. 
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3 The Impenetrable Plaquette Model 



3.1 GENERAL CASE 

The general discretized energy for a flexible phantom triangular fixed- connectivity 



membrane with bending rigidity k is given by [27, 



^ph = \ Yl (V a -rb\-l) + k (l-n a - Hp) , (4) 

<a,b> <",/3> 

where the elastic term is summed over distinct nearest-neighbor pairs of 
monomers a and b and the bending term is summed over all pairs of adjacent 
triangular plaquettes a and (3. In the free energy F ph the parameter e is a 
discrete version of the elastic moduli (the Lame coefficients) of the continuum 
elastic theory pT |, ra 7 is the unit normal to the plaquette 7 and the mean 



lattice spacing has been rescaled to unity without loss of generality. For 
SA membranes F p h must be supplemented by a discretized version of the 
Edwards self-avoidance term introduced in Eq.[l|. 

Without self-avoidance the phantom membrane is crumpled for all tem- 
peratures above the critical crumpling temperature - the microscopic bending 
rigidity is scale-dependent and driven to zero at long wavelength by thermal 
fluctuations. In this regime only the elastic term is relevant in Eq.f|. Under 
the rescaling of the bulk coordinate r —>■ r/ yf]3 we find 

m= \ E (\r a -r b \ 2 - 2y/0\r a -r b \ + const.) , (5) 

1 <a,b> V J 

which reduces to the pure Gaussian model in the high-temperature limit 
(3^0. 

Suppose we can establish that self-avoidance flattens the Gaussian model, 
which corresponds to the elastic term in Eq.|] with vanishing mean lattice 
spacing. Then clearly the SA model Eq.|], together with self-avoidance, will 
also be flat at all temperatures. 

3.2 DISCRETIZATION 

The continuum Edwards model Eq.|l], for membranes {D — 2), is then 

H = J d 2 xd a rd a r + - J d 2 xd 2 x'5 d (f(x) - f(x')) , (6) 
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Figure 1: Two intersecting triangles in three dimensions. In the first case 
the self-avoiding energy is finite, while in the second it is infinite. 



where d is the bulk dimension. 

Let's apply this to piecewise flat surfaces defined by N vertices (f a ) l<a<Nl 
with all vertices, except those at the boundary, being six-coordinated. Any 
point on the surface is within a triangle defined by its three vertices {r a , ff,, r c }, 
and may be parametrized as 

f(a, (3) = ar a + (3r b + (1 - a - 0)r c , (7) 

where a > 0, (3 > and a + (3 < 1. 

The first term in Eq|| is easily discretized as the Gaussian term previously 
discussed 

/ d 2 xd a rd a r = - ^ (f a - r h ) ■ (f a - r b ) , (8) 

1 <a,b> 

whereas for the self-avoiding interaction we have 
J rf 2 xrf 2 x'5 d (f(x) - f(x')) = J d 2 xrf 2 x' f[ <$(r"(x) - r^(x')) (9) 

/d 
da.dpidajdpj J[ 5(r^(x) - r^(x')) , 

where the last sum runs over all distinct pairs of triangles Each term 
here can be explicitly evaluated for any given pair of triangles. Non self- 
intersecting pairs do not contribute to the sum while for self-intersecting 
triangles the result depends on the dimensionality and the precise type of 
self-intersection: 
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• d = 3: Triangles may self-intersect in two different ways, as shown in 
Fig.|l]. In the first case the self-intersection is a line whereas in the sec- 
ond it is a two-dimensional region. One-dimensional self-intersections 
have a finite energy given by 

j dadfidadp 5 {3) (f(a, f3) - f(a, /3)) oc a , (10) 

where s is the intersection length. Two-dimensional self-intersections, 
however, have infinite energies (see Fig.|l]). 

• d — 4: The integral is finite if triangles self-intersect in a point but 
infinite if they intersect in a line or a surface. 

• d > 5 The integral is infinite whenever there is self-intersection. 

The model we treat is an open triangular network containing N vertices 
with free boundary conditions. The energy is 

n = \Y. Va-nW (11) 

<ab> 

provided the triangles defining the surface do not intersect anywhere. If 
two triangles do intersect the configuration is strictly forbidden (its energy is 
infinite). This discretization corresponds to tuning the self-avoiding coupling 
in Eq.[] to infinity (b = oo). It may be worthwhile to consider a weaker version 
of self-avoidance, corresponding to a finite value of b, in which some degree 
of self-intersection is permissible. This will not be treated in this paper since 
the computational cost rises significantly and we do not expect the critical 
behavior to be sensitive to the precise value of b. 

We adopted three distinct geometries, each depicted in Fig.|2|, in our sim- 
ulations. The 3d simulations employed mostly the c geometry, but the other 
two geometries were explored as a test of universality. The 4d simulations 
used the hexagonal geometry a. 

The square geometries b and c, with L vertices on a side, have a total 
number of nodes N = L 2 . The hexagonal geometry a, with L links on each 
of the six sides, has a total number of nodes 

iV = 3L 2 + 3L + l. (12) 
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Figure 2: The three lattice geometries adopted in our membrane simulations, 
where L = 3 for a and L = 5 for b and c. 



3.3 OBSERVABLES and MEASUREMENTS 

A comprehensive understanding of the shape of a fluctuating membrane may 
be obtained from a careful analysis of the shape tensorQ |L9], ffifl , 

^ JV TV iV TV TV 

fc=l 2=1 fc=l fc=l 2=1 

(13) 

where the indices k, I run over the lattice. The trace of S a p is the squared 
radius of gyration R G 

R G = TrS a(3 . (14) 

The growth of the radius of gyration with internal system size L determines 
the Hausdorff dimension du, or equivalently, the size exponent v = 2/dn, 
via Rq ~ L v . Diagonalizing the (symmetric) shape tensor and performing a 
statistical average over independent configurations defines d shape exponents 
(3 a via the finite-size scaling of its eigenvalues X a 

(Xa) ~ Nf" , (15) 

with a — 1, • • • , d. In a crumpled (isotropic) phase all d exponents are equal 
and given by v (i.e. f3 a = v for all a). In a flat phase, on the other hand, 
(3 a = 1 for the two in-plane eigen-directions (say a = 1, 2) and the remaining 

lr The shape tensor is obtained from the standard moment of inertia tensor by subtract- 
ing out the isotropic contribution and is sometimes also mistakenly called the moment of 
inertia tensor in the literature. 
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shape exponents describe the root-mean-square height (h) fluctuations. For 
d = 3, to be specific, /3 3 is the standard roughness exponent ( (< h 2 >~ L 2 ^) 



19, 30 



Further important information on the conformation of the membrane is 
provided by the structure function, defined as 

%"') = ^E<^ ( ^ ) ), (16) 

a,b 

where <f is an arbitrary wavevector. In particular the scaling properties of the 
structure function for wavevectors parallel to the eigenvectors of the shape 
tensor will allow us to determine the key shape exponents of the model. 



4 Simulation Methods 

Simulations of self-avoiding membranes are hampered both by long autocor- 
relation times in updating the embedding coordinates and by the non-locality 
of the self-avoidance constraint. These two factors combined have effectively 
prevented simulations of large enough SA fixed-connectivity membranes for 
a reliable determination of their scaling properties. There has been some 
recent progress though in overcoming both these problems. 

Improved methods for updating the embedding {r(x)} have been applied 
in simulations of a non-SA (phantom) fixed-connectivity membranes. In |31[ 
three different methods were compared: 

• A standard Metropolis updating scheme. 

• Hybrid over-relaxation: make a quadratic approximation to the actionP], 
then apply over-relaxation followed by a Metropolis accept/reject test. 

• Unigrid algorithm: update the membrane recursively on all length 
scales by dividing the lattice into sub-lattices of different sizes and 
apply a Metropolis algorithm to a collective update of those parts. 

The performance of each of these methods - the CPU-cost per indepen- 
dent configuration (Tcpu) - is compared in Fig|| This figure is based on 

2 Note that over-relaxation is exact for vanishing bending rigidity since the action is 
quadratic. 
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Figure 3: The performance of three different methods for updating the em- 
bedding of an isotropic phantom fixed- connectivity membrane. 



the values presented in |y| . For the membrane sizes typically used in Monte 
Carlo simulations, L < 100, both the hybrid over-relaxation and the uni- 
grid algorithm reduced the cost ten-fold compared to a simple Metropolis 
algorithm. Only the unigrid algorithm, however, reduces the dynamical ex- 
ponent z, which measures the volume scaling of the CPU-cost (Tcpu ~ L z ) 
from z = 4 to z ~ 3.8. 

The superior over-relaxation and unigrid algorithms are both easily adapted 
to simulations of SA fixed- connectivity membranes. The real CPU-cost in 
self-avoiding simulations, on the other hand, is the non-locality of the updat- 
ing. Implementing SA as impenetrable plaquettes requires a check verifying 
that the proposed update of the embedding does not lead to intersecting 
triangles. In addition to being non-local, pairwise intersection-checking is 
very time-consuming. It is, however, possible to reduce this CPU-cost by 
a clever implementation of the SA check. We have compared three such 
implement at ions : 

A A comprehensive check; explicitly verifying that no updated triangle 
intersects any other triangle in the lattice by comparing all relevant 
pairs of triangles. 
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Figure 4: The CPU-times for updating a SA fixed-connectivity membrane 
(d = 3) using hybrid over-relaxation with three different methods of checking 
for self-intersection. For each method an estimate of the computational-cost 
exponent is included. 

B Inscribe each triangle in a minimal sphere in the embedding space. This 
has the advantage that pairs of triangles whose spheres do not overlap 
are not compared, and overlapping are quickly identified. 

C Inscribe regions (sub-lattices) of the membrane of different sizes in a 
minimal sphere and apply method B recursively to regions of decreasing 
size. In this way large portions of the membrane are quickly eliminated 
from the checking procedure. 

For simplicity we use these methods in combination with a hybrid over- 
relaxation algorithm as it updates the membrane locally. Although it is not 
a priori clear that the unigrid algorithm, which performs a non-local update, 
should perform any worse than the over-relaxation applied to SA-membranes, 
it would be more complicated to implement the above methods. 

In Fig. 4 we compare the CPU-times for one sweep of hybrid over-relaxation 
for each of these methods. Methods B and C reduce the overall CPU-time 
substantially and, more importantly, both improve the volume scaling of the 
CPU-cost: from t sw ~ L 4A for method A to t sw ~ L 2A for method C. 
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N 


T R G 


IND CONF 


N 


T R G 


IND CONF 


289 


2000 


2500 


2401 


20000 


150 


625 


7000 


1300 


4225 


40000 


80 


1089 


12000 


350 


9025 


71000 


30 



Table 1: Number of independent configurations (IND CONF) for each volume 
for d = 3. 



The simulations were done on several different machines; Pentium II (250 
MHz), IBM SP2 (160 MHz) and a DEC ALPHA workstation cluster. The 
total computational cost of the three-dimensional simulation was equivalent 
to 150,000 SP2 hours while that for the four-dimensional simulation was 
60,000 SP2 hours. The total number of sweeps depends on the volume and 
configurations are stored after a certain fixed number of sweeps. The total 
number of independent configurations for each volume and dimension can be 
read off from Tables [Ij] and |2|. 

We next address the quality of our data. Given an observable O we 
measure, after thermalization, the autocorrelation time To 

ro(T) = 1 - + j:po(t) , (17) 
A t=i 

where the normalized autocorrelation function po(t) is given by 

M-tn.n. /n\2\ 



■ E^r O t O l+t - (O) 



Pott) = — — — 2 ) — . (18) 

and T is a cutoff generally taken to be the largest value for which p does not 
become negative |32| . 



In Tables |1] [d = 3) and ^d = 4) we list autocorrelation times for the 
shape tensor, one of the slowest observables of the membrane. There is a 
striking dependence on volume. 

A cross-check of statistical independence may be performed by studying 
the error bars as a function of bin size. It is well known that the error bar 
for any observable O scales as a function of the bin size as 

a (n b )/a (l) ~ y/(2r + 1) (19) 



3 Note that the autocorrelation times in this table do not give the correct finite-size scal- 
ing exponent z because the larger volumes were simulated with a more effective algorithm 
than the smaller volumes. 
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N MC (xlO 5 ) TH(xlQ 4 ) t Rg IND CONF 



61 


5.5 


5 


8.2(6) 


67073 


127 


8.0 


10 


26.0(9) 


30769 


217 


8.0 


10 


80(10) 


10000 


331 


13.0 


30 


210(35) 


6190 


469 


12.0 


20 


400(40) 


3000 


817 


12.0 


32 


1000(100) 


1200 


1261 


16.0 


40 


3400(700) 


470 



Table 2: Number of Monte Carlo sweeps performed for each volume for 
bulk dimension d = 4, where a Monte Carlo sweep is defined as a combined 
Metropolis and over-relaxation sweep. The number of thermalization sweeps 
is also listed, together with the autocorrelation times of the shape tensor and 
the net number of independent configurations. 

where the bin size n& must be big enough for bins to be statistically indepen- 
dent. To further ensure the statistical independence of our data, we verified 
this dependence in all our measurements. As a concrete example, we plot in 
Fig.[5]the error bar for S22 as a function of the bin size. A clear plateau is ob- 
served corresponding to an error bar a = 0.011. We also find cr(l) = 0.0025. 
^From Eq.|l£] we get r e // = 9.2. Since measurements were taken every 20 
sweeps, we have tr g = 20 and r e // = 184, in good agreement with the result 
given in Table |2|. 

In summary we believe that we have done the requisite checks to ensure 
the statistical independence of our data and we move on to a discussion of 
results. 

5 Measurement of Observables 
5.1 BULK DIMENSION d = 3 

The first observables we analyze are the eigenvalues of the shape tensor Eq.[13|. 
In Fig.^| we plot the distribution of the shape tensor eigenvalues for L = 33 
(N = 1089) and L = 65 (N = 4225). There is a clear first peak, which 
we identify with height fluctuations, and a second peak which we identify 
with size fluctuations. The slight double well structure of the second peak 
is a reflection of the asymmetry of the adopted geometry, as illustrated in 
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_. 200 . 

Bin size 



1000 —. 2000 . 3000 

Bin size 



Figure 5: Plot of a typical error bar as a function of the bin size. The left 
figure magnifies the small size bins of the right one. This particular plot 
correspond to the observable S22 at volume N = 331 (d — 4). Measurements 
were taken every 20 sweeps. 

Fig.0. The overall double peak distribution of eigenvalues corresponds to an 
anisotropic surface resulting from rotational symmetry breaking and is the 
first signal of an orientationally ordered (flat) phase. Indeed, the distribution 
of eigenvalues found here is similar to that found in the flat phase of a 
phantom fixed-connectivity membrane |29[ . 



We next determine the three exponents (3 a from the power law scaling 
given in Eq.|i~5|. The results are shown in Tables ||| and EJ. The column 
corresponding to Fitl is a direct fit to the power law, including all volumes, 
whereas Fit2 excludes volumes N < 625. To investigate the role played 
by the boundary (we have free boundary conditions), we also perform fits 
excluding nodes near the boundary. Excluding nodes beyond 0.9 x L from 
the center, we perform a full fit (Fit3) or a fit removing sizes N < 625 (Fit4). 
We also report similar fits (Fit5 and Fit6) excluding nodes beyond 0.8 x L 
from the center. It is evident that the membrane does not exhibit strong 
edge fluctuations. 

This analysis gives a size exponent v near one and a roughness expo- 
nent ( = 0.65(1). Our results show a slight dependence on the volumes 
included in the fit, indicating that sub- leading corrections are not negligible. 
We repeated the fits using two different parametrizations of the sub-leading 
corrections. This is shown in Table HI. This improved the quality of the fits 
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Figure 6: Distribution of the eigenvalues of the shape tensor for volumes 
N = 1089 and N = 4225 (d = 3). 

considerably. The size exponent v moves towards one and the roughness 
exponent ( decreases slightly. 

Further progress requires study of the structure function, Eq.|TB|. In Fig.^ 
we plot the structure function along the directions of the largest eigenvalue 
for different lattice sizes. For small momentum the structure function is 
monotonically decreasing. For q sufficiently large, however, a series of peaks 
appear. 

^From the small-g region we can extract the critical exponents v and (. 
To do so we plot the structure function along the directions corresponding 





Fitl 


Fit2 


Fit3 


Fit4 


Fit5 


Fit6 


V 


0.945(14) 


0.959(13) 


0.938(13) 


0.948(11) 


0.918(14) 


0.944(13) 


c 


0.655(22) 


0.660(33) 


0.649(19) 


0.655(30) 


0.642(16) 


0.652(23) 



Table 3: Results of the various fits to a simple power law, as described in the 
text, for the size exponent v and the roughness exponent ( (d = 3). 
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fit 


V 


x 2 


c 


x 2 


aNP + b 


0.945(14) 


5.7 


0.655(22) 


0.994 


aNP + blog(N) 


0.959(15) 


5.3 


0.63(1) 


5.9 



Table 4: Comparison of the size and roughness exponents with two distinct 
functional forms of sub-leading corrections (d = 3). The whole lattice is 
included in these fits. 



to the largest and smallest eigenvalues, respectively, as a function of the 
scaling variables qL v and qL^ . Collapse to a single scaling curve is excellent 
for v = 0.95(5) and ( = 0.63(4) as illustrated in Figs.|| {y = 0.95) and | 
(C = 0.63) Q 

The analysis of the structure function for larger values of q is also reveal- 
ing. The picture emerging so far is that a SA fixed- connectivity membrane 
of linear size L is rough at short distances, of the order of L r0 ugh, but globally 
flat. The structure function must exhibit, in this case, peaks at 



thM = TIT + = !,-■■ (20) 

-^rough 

where G\p, are the standard reciprocal vectors for a triangular lattice and 
the factor 7r, as opposed to 2ir, is a consequence of free rather than periodic 



boundary conditions. To test Eq.|20] we plot in Fig.pl] the first (denoted by 
n — 1) and second (denoted by n = 2) peaks of the structure function versus 
1/L. There is clearly a very good fit. Furthermore Eq.|20] implies that the 
ratio of wavevectors g^ =2 /g^ =1 = a/3, which is indeed the case: from the fit 
we find 

bo 46 r- 

2 = ^ = i.77wl.73 = V3, (21) 
b\ 26 

where b n is the slope of the fit (see Fig.|I]). From Eq.[2C] we can also extract 

-trough- 

-trough ~ 8 . (22) 

This result clearly indicates that one must work with membranes of linear 
size L much larger than 8 to effectively eliminate finite-size effects. 

Note that the peaks are damped in intensity with increasing n (see Fig.[7|). 
This may be attributed to the membrane having an effective thickness, as 
given by the roughness exponent. 

4 The quoted errors for these scaling exponents are rather insensitive to the precise 
statistical method by which the quality of the scaling collapse is assessed. 
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Figure 7: Log-Log plot of the structure functions along the direction of the 
maximum eigenvalue of the shape tensor as a function of the wave vector q 
(d = 3). 



The previous analysis convincingly establishes that the membrane is flat. 
Further evidence is provided by visualizing typical snapshots of membranes 
after thermalization. This is shown in Fig.|TT|, where one sees that the mem- 
brane is rough at short distances (I/ roug h ~ 8 by inspection) but flat on large 
scales - certainly "a picture is worth more than thousands of numbers." 
In contrast, a snapshot of a thermalized configuration, with self- avoidance 
switched off (6 = 0), is shown in Fig.|l^. The dramatic effect of self-avoidance 
is striking. As a final check we performed several simulations with a folded 
initial state and observed the subsequent unfolding to the flat phase. 



5.2 BULK DIMENSION d = 4 



We start by examining the shape eigenvalues Eq.|13[ In Fig.[13] we plot the 
distribution of eigenvalues for different volumes. It is clear from the two-peak 
structure that a crumpled phase can be ruled out. The first peak is associated 
with height fluctuations and the second with size fluctuations. These plots 
are in qualitative agreement with those obtained for d = 3. 
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Figure 8: Log-Log plot of the structure function along the direction of the 
maximum eigenvalue of the shape tensor, as a function of the scaling variable, 
for each volume simulated (d = 3). 




Figure 9: Log-Log plot of the structure function along the direction normal 
to the preferred plane of the membrane, as a function of the scaling variable, 
for each volume simulated (d = 3). 
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Figure 10: Fits of the peaks of the structure function to the form q% = 
a + b n /L (d — 3). 

The expectation values of the eigenvalues are plotted as a function of 
volume in Fig.[L4|. Both Ai and A 2 increase rapidly with volume, while A 3 and 
A 4 grow slowly. 

To compute the exponents /3 a , we need to fit the data to the functional 
form in Eq.[15|, which has three free parameters. 

Finite size effects are significant for the range of volumes we analyzed 
(N = 127 to 1261). To evaluate the importance of sub-leading corrections 
to scaling, we once again fit the data to two distinct functional forms. The 
results are shown in Tables [| and |6|. The size exponent spans the range 

v = 0.82(4)- 0.90(1) , (23) 

Although sub-leading corrections significantly affect the extracted scaling 
behavior, it is clear that the size exponent v is tending to 1. To substantiate 
this result requires extensive simulations of larger volumes. 

In a flat phase, the scaling of the third and fourth eigenvalue is associated 
with the roughness exponent. The roughness exponent falls in the range 

C = 0.65(1) - 0.79(1) , (24) 

where again, the errors ignore systematic effects. The roughness exponent 
is large, so it implies that the membrane, despite being flat (or almost flat) 
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Figure 11: A snapshot of a thermalized self-avoiding configuration for volume 
N = 9025. The membrane is flat over long length scales but rough on short 
scales (d — 3). 




Figure 12: A snapshot of a thermalized configuration (d = 3) for a phantom 
membrane (b = 0), with vanishing bending rigidity, for comparison with 
Fig.|TT] (N = 4225). 
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Figure 13: Distribution of eigenvalues for three different volumes N = 331, 
N = 817 and N = 1261 (d=4). 



overall, is very rough. For completeness, we quote the roughness exponent 
associated with the smallest eigenvalue; it is ( u = 0.62(7) — 0.74(2). It is 
somewhat smaller than Eq.^ but compatible with it. Whether these two 
exponents are the same is something we cannot establish, and larger volume 
sizes are clearly called for. 

The other observable we use to study the properties of the membrane 
is the structure function, defined in Eq.|i~6|. As noted earlier, the scaling of 
the structure function along the direction defined by the largest eigenvalue is 
directly related to the size exponent. The strategy is to match the structure 
function as a function of the scaled variable qL u , where v is the size exponent. 
The result of the matching is 

v = 0.9(1) , (25) 
where the error bar needs further explanation. The matching of the structure 
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Figure 14: Plot of the eigenvalues as a function of the volume (d=4). 



function for small values of q is acceptable for values v > 0.8. Near v = 0.8 
the range of q for which matching is found is somewhat small. This range 
increases with u, but if v is chosen too large (say, v > 1), the small-g matching 
deteriorates and is poor. We note that although the matching for v = 1 is 
still very good, the best overall value for v is v = 0.9. The error bar quoted 
in Eq. |25] is a conservative estimate embracing the range of acceptable v. 

We also computed evaluated the structure function along the directions of 
the third and fourth eigenvalue of the shape tensor. This yields a roughness 
exponent (see Fig.|i~6|) 

C = 0.60(7) , (26) 

in good agreement with the estimates given in Tables [S] and []. The error bar 
is a rough estimate indicating the range of exponents that yield acceptable 
scaling. 

The situation for the plaquette model for d > 4 is much easier to analyze. 
Since a two-dimensional surface almost never self-intersects in bulk dimen- 
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aN^ + b 


P 


a 


b 


x 2 


Ai 


0.833(8) 


0.0172(9) 


0.47(2) 


0.24 


A 2 


0.81(2) 


0.013(2) 


0.20(3) 


0.61 


A 3 


0.65(1) 


0.0119(9) 


0.166(7) 


0.22 


A 4 


0.620(7) 


0.0093(6) 


0.112(5) 


0.10 



Table 5: Fits to finite-size scaling of the form aN 13 + b for the eigenvalues of 
the shape tensor (d — 4). 



aN^ + blog(N) 


P 


a 


b 


x 2 


Ai 


0.90(1) 


0.010(1) 


0.139(6) 


0.48 


A 2 


0.85(2) 


0.0091(2) 


0.060(7) 


0.60 


A 3 


0.788(6) 


0.0038(1) 


0.0574(1) 


0.10 


A 4 


0.74(2) 


0.0031(4) 


0.038(1) 


0.22 



Table 6: Fits to finite-size scaling of the form aN 13 + b log(iV) for the eigen- 
values of the shape tensor for (d=4). 

sions five and above, self- avoidance is irrelevant and the membrane is always 
crumpled. The radius of gyration for any d > 5 is given by 

R 2 G ~ log(TV) , (27) 
which corresponds to v — 0, the Gaussian result. 

6 Summary and Conclusions 
6.1 SUMMARY 

In this section we summarize the most important results obtained from the 
detailed analyses presented earlier. The impenetrable plaquette model for 
d — 3 is flat for all temperatures with critical exponents 

v = 0.97(4) C = 0.63(3) , (28) 

with the error reflecting the range of estimates obtained by the different 
methods described. The membrane is a flat object at large distances but a 
very rough one for characteristic linear sizes L ~ 8. The values for the critical 
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Figure 15: Log-Log plot of the structure function (d=4) along the direction 
of the maximum eigenvalue of the shape tensor for each volume simulated as 
a function of the scaling variable qL u (with v = 0.9). 



exponents are in agreement with the results from BS models (see @ for a 
review). The roughness exponent of the flat phase of phantom membranes 
is given by p| 

C = 0.64(2) , (29) 

which agrees very well with the result obtained above. We conclude that there 
is a single flat phase for fixed-connectivity membranes, describing either a 
phantom membrane at large bending rigidity or a self-avoiding membrane at 
any non-negative bending rigidity. 

Embedding dimension d = 4 has also been studied. Here we find very 
strong evidence that, for all temperatures, there is once again only a flat 
phase, with exponents 



i/ = 0.9(1) C = 0.65(10). 



(30) 



In this case the quantitative values for the exponent are not as accurate as 
we have for d = 3 since the lattice volumes considered were not as large 
and finite size effects give rise to systematic errors. It would be desirable to 
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Figure 16: Log-Log plot of the structure function (d=4) along the direction 
of the minimum eigenvalue of the shape tensor as a function of the scaling 
variable qL/* (with ( = 0.6). 



simulate larger volumes to narrow down the exponent v ~ 1. Gathering all 
the evidence acquired from the analyses carried out in this paper, however, 
we find it very unlikely that the size exponent is not one. 

6.2 CONCLUSIONS 

Our study establishes that the phase diagram of fixed-connectivity mem- 
branes is very simple and consists only of a flat phase for both d = 3 and 
d — 4, and crumpled phases only for d > 4. In order to rigorously compare 
these theoretical results with experiments, one should consider the effect 
of topological defects | 55fl , which play an important role in any crystalline 



phase. One would not expect defects to alter the long-distance properties 
of the model, however, within the range of temperatures for which the crys- 
talline phase prevails, since the overall integrity (unbroken connectivity) of 
the mesh combined with self-avoidance are the key triggers of the flat phase. 
These properties are retained even in the presence of topological defects. This 
overall picture is consistent with the current experimental situation - no re- 



26 



alization of a fixed-connectivity membrane has ever been seen in a crumpled 
state. The observed flat phases have exponents consistent with the numerical 
results of this paper. 

Analytical results for SA fixed-connectivity membranes are reviewed in 
H. Here we simply note that the e-expansion, at zero bending rigidity, 
predicts a unique infra-red stable fixed point and no phase transitions, in 
agreement with our simulations. The current value of the size exponent 



within the e-expansion, computed to two- loop order in [34, p5j (see B for a 



review), however, definitely seems to be less than one. It would be of great 
interest to know if the inclusion of higher-orders in the expansion pushes u 
towards the flat value of one. 

It is also important to remark that the previous results assume repulsive 
potentials among the monomers. When this restriction is dropped and at- 
tractive forces are considered [T9| , the picture actually changes, and compact 



phases |36| seem to be possible (see also |37|, |38| . This very interesting possi- 
bility may naturally occur in some systems, so it will certainly be the subject 
of subsequent work. 

Our results also have implications for other models. It is well known that 
anisotropic fixed-connectivity membranes possess an intermediate tubular 



phase |39| |40[ . Since the tubular phase is, very roughly speaking, crumpled 
in one dimension only, it has been argued that it may survive the incorpora- 
tion of self-avoidance HI, 42], |43| . Anisotropy cannot be readily introduced 



in our model, since if we tune the self-avoidance coupling to zero, we au- 
tomatically sit in the crumpled phase where anisotropy is known to be an 
irrelevant perturbation. Since adding anisotropic extrinsic curvature to the 
self-avoiding membrane will only flatten the membrane even more, this can- 
not produce a tubular phase either. One should therefore consider different 
discretizations for the self-avoiding case than the ones that are suitable for 
the phantom model discussed in [fi(] . 

The understanding of the physical properties of fixed- connectivity mem- 
branes are also very important in constructing realistic models of full-fledged 
cell membranes. The next step in this very exciting goal would be to de- 
scribe a model of a coupled fluid and fixed- connectivity membrane. Since 
one can safely assume (provided no attractive forces are present) that the 
fixed-connectivity membrane is flat, this property alone significantly con- 
strains the different effective theories that need to be considered. This is just 
one of the many exciting problems that physicists will need to tackle in the 
near future. 
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